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Abstract. We develop a Coulomb gas description of the critical fluctuations in 
the fully packed loop model on the honeycomb lattice. We identify the complete 
operator spectrum of this model in terms of electric and magnetic tjecior-charges, 
and we calculate the scaling dimensions of these operators exactly. We also study 
the geometrical properties of loops in this model, and we derive exact results for 
the fractal dimension and the loop size distribution function. A review of the many 
different representations of this model that have recently appeared in the literature, 
is given. 
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1. Introduction 

Loop models, loosely speaking, are statistical models which have as basic building 
blocks loops that run along the bonds of a two-dimensional lattice. In order to define 
a loop model one assigns Boltzmann weights to the different loop configurations. This 
is usually implemented by assigning weights to the different vertex configurations 
allowed by the loop model, and a weight to the loop as a whole; this loop-weight is 
usually referred to as the fugacity. 

Loop models have attracted attention recently as representations of certain 
exactly-solvable vertex models which can be used to construct restricted solid-on-solid 
models Q, some of which admit an off-critical extension ||. They are also particularly 
simple examples of models that allow for Monte-Carlo simulations with non-local 
loop updates, which have been recently studied as algorithms that reduce critical 
slowing down || . In a completely different setting, loop models appear as space-time 
diagrams in the path integral representation of one-dimensional quantum spin chains, 
and many quantities defined in terms of the spins can be reexpressed in the language 
of loops (J] . For instance, it can be shown that the spin-spin correlation function in the 
antifcrromagnetic Heisenberg spin chain, is given by the loop correlation function in 
the appropriate loop model Q . The loop correlation function measures the probability 
that two points on the lattice belong to the same loop. 

In two-dimensional classical spin models loops are typically encountered as domain 
boundaries, e.g., Bloch walls in the Ising model, or as graphical representations of 
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high-temperature expansions. Recently, Cardy |5| has calculated different geometrical 
properties of cluster boundaries in the 0(n) model on the honeycomb lattice, using 
the loop representation of this model. 

In this paper we study the fully packed loop (FPL) model on the honeycomb 
lattice. The FPL model was introduced by Reshetikhin Q , and independently by Blote 
and Nienhuis 0, as the zero-temperature limit of the 0(n) model. Using numerical 
transfer-matrix methods, Blote and Nienhuis were able to to show that this loop model 
defines a universality class distinct from the previously studied low temperature phase 
of the 0{n). Exact values of the critical exponents and the conformal charge of the 
FPL model were subsequently determined by Batchelor et al. ||, who found a Bethe 
ansatz solution of the model. 

In the FPL model nonintersecting loops are placed along the bonds of a honeycomb 
lattice so that each vertex of the lattice is covered by a loop. The partition function 
is given by 

Z FPL =J2n N ^ (1) 

g 

where N is the number of loops in the fully packed configuration Q, and n is the loop 
fugacity. 

This model undergoes a phase transition as a function of the loop fugacity n. 
For values of n approaching zero, configurations with a small number of big loops 
are favored; in the limit n — ► a single loop covers the whole lattice.^ Loops of all 
sizes will be present on the lattice as we increase n. This is equivalent to having a 
diverging correlation length in the system ||, and the model is critical with power 
law correlations. At n = 2 the FPL model undergoes a Kosterlitz-Thouless type of 
transition |j, [l0| into a long-range ordered state, in which there exists a largest loop 
on the lattice which is roughly the size of the correlation length. In the n — >• oo limit 
the fully packed loop model is perfectly ordered, with all the loops having the minimal 
length of six, and occupying one of the three sublattices of hexagonal plaquettes. 

Here we focus our attention on the FPL model along the critical line (0 < n < 2), 
which was also the focus of the above mentioned numerical transfer-matrix study, 
and of the Bethe ansatz solution. Using a nested Bethe ansatz Batchelor et al. || 
calculated the scaling dimensions of the "watermelon" operators along the critical 
line. The watermelon correlation function is defined as the probability that m loop 
segments meet in the neighborhood of two points separated by r Here we 

rederive the same results from a Coulomb gas approach in which the loop model 
is mapped to an interface model. In the interface representation "watermelon" scaling 
dimensions become associated with vortices whose topological charges are vectors in 
the triangular lattice. Furthermore, this approach allows us to identify the complete 
operator spectrum of the FPL model and make contact with known results from 
conformal field theory. In particular, we calculate exactly the temperature dimension 
found numerically by Blote and Nienhuis |7) , that does not appear in the Bethe Ansatz 
solution of Batchelor et al. We also identify defect configurations in the loop 

model that generalize the "watermelon" configurations, and we calculate the critical 
exponents associated with them. 

The main shortcoming of the Coulomb gas approach, in general, is that it usually 
relies on some exact information about the model which can then be used to calculate 

f This limit was studied by Batchelor et al. [^| who obtained exact results for Hamiltonian walks on 
the honeycomb lattice. 
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the value of the renormalized coupling [ p2[ . Once the coupling is known all the 
exponents can be calculated exactly. We will show that in the FPL model the coupling 
can be determined exactly by identifying the marginal operators in this model. The 
existence of these operators is required by certain consistency conditions placed on 
the conformal field theory which describes the scaling limit of a lattice model; this 
was discussed at length by Dotsenko and Fateev [jl3| . 

Many of our results for the critical exponents of the FPL model have been found 
previously from a Bethe Ansatz solution [gj. Our main motivation for pursuing the 
Coulomb gas approach is its relative simplicity, and the geometrical interpretation of 
the operator spectrum of the FPL model, which it offers. Furthermore, this approach 
allows us to identify the conformal field theory that describes the scaling limit of the 
FPL model, which can then be used to study the critical properties of this model in 
detail, using the many conformal techniques at our disposal. 

This paper is organized as follows: in section ^ we have collected the different 
known representations of the FPL model, and mappings from one to the other are 
made explicit. In section || we introduce an effective field theory of the FPL model, 
and section [| is devoted to the construction of the associated Coulomb gas and the 
calculation of the scaling dimensions of different operators. These results are used in 
section || to calculate the geometrical exponents for loops in the FPL model. 



2. Representations of the FPL model 

The FPL model has many different representations, some of which have been 
independently studied. Here we review the mappings between the different 
representations . 

The FPL model for < n < 2 is equivalent to the three-colouring model on the 
honeycomb lattice introduced by Baxter [ f[4| . The three-colouring model is defined 
by colouring the bonds of the honeycomb lattice with three different colours, say A, 
B, and C, in such a way that no two bonds of equal colour meet at a vertex. [f] For 
n = 2 each colouring is given equal statistical weight. If we choose any two colours, 
say B and C, then the bonds coloured with these two colours form a fully packed loop 
configuration on the honeycomb lattice. Each loop can be coloured with alternating 
colours B and C in two ways (B-C-B ... or C-B-C . . .), and is therefore assigned a 
fugacity n — 2. The FPL model away from the n — 2 point can also be mapped 
to a colouring model, but now the weights of the different colourings will have to be 



modified; see section 4.2 



If we consider the three colours as Potts spins placed at the centers of the bonds 
of the honeycomb lattice, then the three-colouring model describes the ground state 
of the three-state antiferromagnetic Potts model on the Kagome lattice studied by 
Huse and Rutenberg |jl5). The Potts model is defined by the Hamiltonian (energy 
functional) 

H = \J\J2 S *w > ( 2 ) 

<ij> 

where the sum goes over nearest-neighbors, and the spins <Ji live on the vertices of the 
Kagome lattice. At zero temperature the only allowed states are ones where on every 
triangular plaquette all three spins are present. This ground state manifold is critical 



f These type of graph colourings are known in the mathematics literature as edge colourings. 
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Figure 1. One of six symmetry related ideal states of the three-colouring model. In 
an ideal state all the plaquettes are coloured with two colours only. The microscopic 
heights z are defined at the centers of the plaquettes, and the change in z, when 
going from one plaquette to the neighboring one, is determined by the colour of the 
crossed bond. The ideal state is macroscopically flat, in the sense that the variance 
of the microscopic height is minimal. 



in the sense that correlation functions of the spins decay with distance as power laws 

The three colouring model can be mapped to a solid-on-solid model that describes 
a two dimensional interface in four spatial dimensions |H| . This is accomplished 
by placing a two component microscopic height z at the center of each plaquette of 
the honeycomb lattice; see figure [l| The change in z when going from a plaquette to 
the neighboring one is given by A, B, or C, depending on the colour of the bond that 
is crossed; the vectors A, B, and C point to the vertices of an equilateral triangle, 

where we have chosen the normalization for later convenience. We also adopt the 
convention that the microscopic height increases by A, B, or C when going clockwise 
around an up pointing triangle of the dual lattice (i.e., around a vertex of the 
honeycomb lattice in the shape of the letter Y); see figure [I]. Up to an arbitrary 
choice of a single height, say at the origin, the mapping of the colouring to the heights 
is one-to-one. Each allowed height configuration is given equal statistical weight. 

It has been recently shown by Di Francesco and Guitter |l(J that the three 
colouring model can be mapped to a folding model of the triangular lattice, introduced 
by Kantor and Jaric ]l7| . The allowed configurations in the folding model are given 
by all the possible complete foldings of the triangular lattice.^ A folding configuration 
can be specified by giving the direction (+ or — ) of the normal to each elementary 
triangle in the folded state. Now if we place the three colours, A, B, and C, on the 
bonds of every elementary triangle in a clockwise (+) or anticlockwise (— ) fashion 
we obtain a three colouring configuration of the dual honeycomb lattice. This is a 
six-to-one mapping, since for a given folded configuration one is free to choose one of 

f This model is a special case of a more general folding model due to Shender et al. Jllj], which is 
equivalent to the ground states of the antiferromagnetic Heisenberg model on the Kagome lattice. 
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six colour configurations around a single triangle, which then fixes all the rest. 
3. Effective-field theory 

In this section we propose an effective field theory for the long wavelength fluctuations 
of the interface model, which is one of the representations of the FPL model discussed 
in the previous section. Here we focus on the n = 2 case which is equivalent to the 
three-colouring model with equal statistical weight for all the colourings, and extend 
to n < 2 in the following section. 

We motivate the long- wavelength theory of the interface model by a coarse-graining 
procedure of the microscopic heights z, which is implemented as follows:|| First, we 
define the ideal states which we use to coarse-grain the three-colouring model. Ideal 
states are edge colouring states in which every elementary plaquette of the honeycomb 
lattice is coloured with two colours only, figure |l|. These states are flat, in the sense 
that they have the smallest variance of the microscopic height, and we argue that the 
free energy of the colouring model || is dominated by fluctuations around the ideal 
states. Namely, the smallest change on the lattice, that is allowed by the constraints 
of the three-colouring model, is an exchange of colours along a loop of alternating 
colour (eg. C-B-C ... to B-C-B . . .). The ideal states maximize the number of loops 
that allow for these loop exchanges, and it is this property that selects them out. This 
entropic selection effect is close in spirit to the "order by disorder" effect introduced 
by Villain |2G[ ]. In the FPL model ideal states are the ones selected in the n — + oo 
limit, while in the folding model these are states where the triangular lattice has been 
folded to a single triangle. 

Second, we divide the honeycomb lattice into domains so that each domain 
represents a different ideal (flat) state. To each domain we assign a coarse grained 
height h, which is equal to the microscopic height averaged over the domain; h = (z). 

The coarse grained heights associated with the six different ideal states (one for 
every permutation of the three colours) form a honeycomb lattice which we call the 
ideal state graph X; see figure |^. The side of the elementary hexagon of X is 1/3, in 
the units chosen for the vectors representing the colours, equation (||). Nodes of X 
that correspond to the same ideal state form a triangular lattice with an elementary 
triangle of side 1. This lattice we call the repeat lattice TZ; in the following section 
we will show that vectors b € 1Z are the magnetic vector-charges in the Coulomb gas 
associated with the FPL model. 

Finally, we consider the continuum limit of the interface model in which the discrete 
heights, defined over different ideal state domains, are replaced with a continuously 
varying height field h(r) = (hi(r), h 2 (r)). The dimcnsionless free energy / of the 
interface, which is entropic in origin, is assumed to be of the form 

/ = J dV [^(iV/ixl 2 + |V/i 2 | 2 ) + V(h)] , (4) 

where 2ng is usually referred to as the stiffness. V(h) is a periodic potential with the 
periodicity given by the ideal state graph, i.e., 

V(h + X) = V(h) . (5) 

f Details of the height construction for the n = 2 FPL model, as well as for other critical ground 
states, can be found in 

f The free energy of the three-colouring model is purely entropic in origin, in the sense that the 
partition function is simply equal to the number of different edge colourings. 
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Figure 2. The ideal state graph of the three colouring model is a honeycomb 
lattice in height space: each vertex is associated with a particular ideal state, and 
the six different ideal states form a hexagonal plaquette. The ideal states are labeled 
by the colour configuration (crj , CT2 , CT3) of the bonds around a common vertex of the 
coloured lattice. The vertices in the ideal state graph that correspond to the same 
ideal state (say (C, B, A)) form a triangular lattice which is the repeat lattice of the 
three-colouring model. 



The free energy / defines an effective field theory of the three-colouring model; the 
assumption being made is that it correctly describes the long-wavelength fluctuations 
of the microscopic height z. The periodic potential V(h), which is usually referred to 
as the locking potential, favors the heights to take their values on X, while the first 
term is the contribution to the free energy due to fluctuations around the flat ideal 
states. Therefore, the assumption that the effective field theory of the three-colouring 
model is given by equation (Q) is directly related to the intuitive idea put forward 
earlier, that the free energy of the three-colouring model is dominated by fluctuations 
around the ideal states. 

The locking potential is periodic with the periodicity of T. Thus, the three- 
colouring model, in its interface representation, undergoes a roughening transition 
for some value of the coupling g — g r If the coupling g satisfies g < g r , then 

the locking potential in equation (Q) becomes irrelevant, in the renormalization group 
sense, and the three-colouring model is described by a Gaussian model with a free 
energy 

f = ng /dVdV^p + lV^I 2 ) . (6) 



In the case that the locking potential is relevant (g > g r ), the three-colouring model 
will lock into long range order in one of the ideal states. This would imply a finite 
correlation length in the FPL model, roughly the size of the largest loop in the system. 
We will see in the following section that for the the three-colouring model which is 
equivalent to the n = 2 FPL model, g is equal to g r , so the interface is at the roughening 
transition (i.e., the locking potential is marginal). For values of the fugacity n < 2 it 
will be shown that g < g r , but due to the presence of the background charge the the 
locking potential remains marginaljf] 

f Marginal roughness seems to be a general property of critical loop models in two dimensions p^ ] . 
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4. Coulomb Gas 

Here we develop the Coulomb gas description of the FPL model based on its interface 
representation. We determine the spectrum of possible electric and magnetic charges 
and we calculate the scaling dimensions of operators associated with them. These are 
compared to recent numerical results Q| , and to results from a Bethe Ansatz solution 
of the FPL model || . We first consider the n = 2 FPL model which is described by a 
simple Gaussian field theory, equation (|^) . The n < 2 case is treated by perturbing the 
n — 2 theory with an integrably marginal operator p3[ , and introducing a background 
charge in the Coulomb gas . 

4-.1. FPL model at n = 2 

In constructing the effective field theory of the three-colouring model, which maps 
to the n = 2 FPL model, we found that the height field h is compactified on the 
triangular lattice TZ. Therefore, any local lattice operator $(r ) uniform in the ideal 
states is periodic in height space, and it can be written as a Fourier series 

*(r ) = $ G e i2 " 5 K(F) . (7) 

Here TZ* is the lattice dual to TZ, i.e., G ■ b £ Z for any two vector-charges b £ 1Z and 
G £ TZ* . The Gaussian field theory, equation (Q), describes the vacuum phase of a 
two-dimensional Coulomb gas of electric (G) and magnetic (6) vector-charges [Q . 

The scaling dimension of $(r ) is equal to the scaling dimension of the most relevant 
vertex operator exp (i27rG • h(r)) appearing in its Fourier expansion. The scaling 
dimension x(G), of a vertex operator, can be easily determined from its two-point 
correlation function, 



e 



i27rG-/i(r) e -i2?rG-/i(0)^ ^ r ~2x{G) ^ 



Namely, using a general property of a Gaussian distributed random field 

,^d-h(¥) e -i2nG-h(0)\ = p -i([27rG-(K(r)-h(0))] 2 ) 



e -§U3f«-W'v--'»w;r; ) (9) 



and the height-height correlation function calculated in the Gaussian field theory 
defined by equation (^|) (for r»a; a is the lattice spacing) 

((h t (r) - ^(O)) 2 ) = const + ^- In \r\ , (10) 

we find: 

x{G) = f |G| 2 . (11) 

Operators with a non-zero magnetic charge b £ TZ can be associated with a vortex 
configuration of the height field, or a violation of the edge colouring constraint in the 
three-colouring model; see figure |^. The topological charge of the vortex is given by 
6, i.e., the height mismatch around a vertex of the honeycomb lattice, and it can be 
calculated using the height rule introduced in section ||. 




Figure 3. Elementary defects in the three-colouring model are associated with 
loops of alternating colour; exchanging the two colours (B and C) along one half of 
the loop (shown in bold) will generate defects (circled), that is violations of the edge 
colouring constraint, at the two ends. In the interface representation these defects 
become vortex-antivortex configurations of the height. 



For a magnetic- type operator the scaling dimension x(b) follows from tl 
expression for the (dimensionless) interaction energy of a vortex-antivortex pair jl 
(for r>«) 

= g\b\ 2 In r + const . (12) 



E 



Here ±6 are the topological charges of the vortices and r is the distance between 
them: the above expression for E- mt follows from the Gaussian form of the free energy, 
equation (^), and is simply the two-dimensional version of Coulomb's law. The vortex- 
antivortex correlation function is given by the Boltzmann factor 



1 



„2x(b) 



where x(b) is fixed by the coupling g, 

X (b) = 1 16| 2 . 



(13) 



(14) 



Equations (11) and ( |14| ) specify the complete spectrum of scaling dimensions in 
the n = 2 FPL model. Next we turn to the problem of identifying specific correlation 
functions in this model associated with these exponents. 

The magnetic dimension Xb., introduced by Blote and Nienhuis 0, governs the 
probability that two points on the honeycomb lattice lie on the same loop of the FPL 
model. The loops in the n = 2 FPL model are contour loops of a particular component 
of the height field. This is an important observation which will lead us to the exact 
value of the coupling g, and we clarify it further. 
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Say we choose the BC loops in the three colouring model to represent the loops in 
the FPL model. Now consider the points at the centers of the hexagonal plaquettes 
along one side of a BC loop. These points are separated by A coloured bonds, and 
consequently the component of the microscopic height z orthogonal to A is unchanged 
as we go along the loop. Recently it has been argued that the exponent associated with 
the correlation function that measures the probability that two points belong to the 
same contour loop of a random Gaussian surface is independent of the stiffness {2ng), 
and equal to 1/2 j24|. Therefore, we conclude, the magnetic dimension in the n = 2 
FPL model is Xh = 1/2. The numerical result of Blote and Nienhuis is x^ = 0.470(1) 
(see Table I in 0), and x^ — 1/2 was also found from the Bethe Ansatz solution ||. 

In the interface representation of the three-colouring model the magnetic dimension 
can also be associated with a vortex- antivortex correlation function with the magnetic 
vector-charge 

b h =B-C = (0, 1) . (15) 

This comes about in the following way. The correlation function that measures the 
probability that two points separated by r belong to the same BC loop is Z(r)/Z . 
The restricted partition function Z(f) is simply the number of colourings with a BC 
loop passing through and f, while Zq is the total number of colourings. Now, if we 
exchange the two colours on the loop along one half of the loop going from to r, then 
this will generate a vortex and an antivortex with charges ±&h, &t these two points. 
The magnetic dimension Xh is therefore 

Xh = x(b h ) . (16) 

Earlier we found x^ = 1/2, and from equation ( |l4| ) we can now calculate the exact 
value of the coupling in the n = 2 FPL model: 

.9 = 1- (17) 

Huse and Rutenberg [[TJ) deduced the value of the coupling constant g for n = 2 
from the exact solution of the three-colouring model due to Baxter ll4|] . They showed 
that from Baxter's solution one can conclude that the interface model is exactly at 
its roughening transition. In standard scaling terms this means that the dimension of 
the locking potential, V(h) in equation (0), is 2. This on the other hand leads to the 
equation x(Gy) = 2, where Gy is the smallest electric vector-charge appearing in the 
Fourier expansion (ffl)of the operator V(h); Gy is also the second shortest vector in 
1Z* . From equation ( fil| ) and \Gy\ = 2 the value of the coupling, g = 1, follows. 

Another dimension that was measured by Blote and Nienhuis is the temperature 
dimension which is associated with a vertex of the honeycomb lattice not covered by 
a loop. An uncovered vertex becomes in the three-colouring model a defect with the 
same colour on all three of its surrounding bonds JToJ • in the interface representation 
this becomes a vortex with a topological chargefj] 

b t = 3A, (18) 

which is the second shortest vector in 1Z. Hence, using equation (|lj) for the dimension 
of a magnetic-type operator, and equation (|l7]) for the value of the coupling, the 
temperature dimension is 

xt = x(b t ) = 3/2 . (19) 

f Here and throughout we take the loops in the FPL model to be BC coloured. Of course any other 
choice of colours would lead to the same results. 
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For comparison, the numerical transfer matrix result is x t = 1.46(1) (see Table I in 
0), and it is in good agreement with our exact result when systematic errors are 
taken into account j?J. This dimension was not calculated by Batchelor et al. ||. 

One possible generalization of the loop correlation function is the so called 
watermelon correlation function which measures the probability of having m loop 
segments meeting at two point separated by r flTl[ ; m = 2 would then be the loop 
correlation function. The central result of the Bethe Ansatz solution of the FPL model 
by Batchelor et al. ^] was the calculation of the watermelon dimensions x m . Here we 
show how they can be calculated from the Coulomb gas. 

In a way that is completely analogous to the above analysis of the loop correlation 
function, the watermelon correlator becomes a vortex-antivortex correlation function. 
Therefore, in order to calculate x m we need to determine the appropriate magnetic 
vector-charges b m . We divide this calculation up into four steps: 

(i) m = 1 — This corresponds to having a single BC loop segment between points 
and r, which in turn implies that the colour configuration at is {A, A, B}. 
Using the height rule around we find 

b 1 = 2A + B = (^-,±). (20) 

(ii) m = 2 - In this case there is a BC loop originating at and ending at r. The 
associated topological charge was found earlier, equation (|l5|), 

b 2 = b h = (0,1) (21) 

(iii) m = 2k - Here we have k BC loops originating in the neighborhood of 
and ending in the neighborhood of r. The total topological charge around the 
endpoints is 

b 2k = kb 2 • (22) 

(iv) m = 2k — 1 - This case can be thought of as having k — 1 BC loops and 
an additional BC segment, originate in the neighborhood of 0, and end in the 
neighborhood of r. The total topological charge is 

b 2k -i = {k - l)b 2 + h , (23) 

where b\ and b 2 are given by (^0|) and (pl|), respectively. 

Using the exact value of the coupling, g = 1, and the above calculated magnetic 
vector-charges, from equation ( |l4| ) we find for the watermelon dimensions: 

x 2 k = -jk 

x 2k -i = i(fc 2 -fc + l) . (24) 

These exponents form a special case of the more general spectrum of exponents 
associated with all the possible magnetic vector charges in the repeat lattice foj.fc = 
j{A-B) + k{A-C) eTZ 

Xj, k = x(b jtk ) = - (f + k 2 + jk) , (25) 
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where we have once again made use of equation (|1J) . 

These charges are associated with defects in the FPL model, in which the full- 
packing constraint is violated locally. Each defect can be contained in a non- 
selfintersecting polygon on the triangular lattice, dual to the original honeycomb 
lattice. The edges of the honeycomb lattice cut by the polygon will be called defect 
edges. 

The topological charge of a defect in the solid-on-solid version of the model is the 
vector sum of the height differences measured along the polygon in, say the clockwise 
direction. This trivial procedure yields via equation ( p5| ) the exponent governing the 
spatial decay of correlations between two oppositely charged defects. In the three- 
colouring model the charge of a defect is simply that of the corresponding height 
configuration. 

In the loop version of the model two-defect correlations are defined by the 
requirement that specified defect edges of the two defects be connected by a loop 
segment, and that the remaining edges be empty. Like everywhere else in the lattice 
the empty edges correspond to the colour A and the occupied ones to B or C. The 
possibility of a loop segment connecting defect edges of the same defect must be 
excluded. This can be done by maximizing the charge of the defect using the choice 
between B and C . Thus the corresponding topological charge of the one defect can be 
found by associating the vector A or —A to the empty defect edges and B or — C to 
the occupied ones. For the other defect the same rule applies with B and —C replaced 
by C and —B. 

It may be noted that the two sublattices of the honeycomb lattice are not equivalent 
when positioning a defect. For instance according to the above rule the charge of a 
vertex with three empty edges is 3A on the one sublattice but — 3A on the other. Two 
such defects have a total charge zero only if they are placed on different sublattices. 
With periodic boundary conditions correlations functions between two defects that do 
not have opposite charge are zero. However, with open boundary conditions there are 
non-zero two-point correlations between different magnetic defects; for defect charges 
b and b' the exponent is x(b) + x(b) — x(b + b'). 

4.2. FPL model for n <2 

Here we extend the Coulomb gas description of the FPL model for loop fugacities 
n < 2. We show that this can be accomplished by introducing a background vector- 
charge. The effect of the background charge is a lowering of the conformal charge, 
and a shift in the scaling dimensions found above. 

Let us first introduce, in the three-colouring model, the staggered chirality x(f)- 
This operator takes two values: it is +1 (—1) if the colours go clockwise around the 
vertex r on the even (odd) sublattice of the honeycomb lattice, and —1 (+1) if the 
colours go anti-clockwise. For n < 2 the free energy (||) has an extra imaginary bulk 
term conjugate to the xffl'i this term was introduced in the colouring model by Baxter 
Q . The effect of the bulk term iA J d 2 rx(r), is to assign a phase factor exp(±iA) 
every time a loop in the FPL model makes a left or a right turn. This will have the 
effect of assigning to each loop a weight 

n = e~ i6A + c i6A = 2 cos(6A) , (26) 

due to the fact that the difference between the number of left and right turns, when 
walking along a closed loop on the honeycomb lattice, is six. By inspection of the ideal 
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state graph in figure Q we conclude that the most relevant vertex operator appearing 
in the Fourier expansion of x{^) nas an electric vector-charge G x , which is the second 
largest vector in 1Z* (same as for the locking potential). Therefore at fugacity n = 2 
the scaling dimension of the chirality is 2, i.e. it is marginal. We expect that its 
only effect on the Gaussian action / in equation (^) will be an isotropic change of the 
coupling g. | 

Other then the marginal term, the shift of the loop fugacity away from n = 2 will 
also generate a term in the free energy (Euclidean action) which couples the height 
field to the scalar curvature j25j. This can be most readily understood by taking the 
FPL model to be defined on a cylinder of circumference L. Namely, a seam running 
along the length of the cylinder has to be introduced in order to give the correct 
fugacity (n) for loops winding around the cylinder. A bond of colour a S {A,B,C} 
which crosses the seam gets an extra factor exp(i27r£?o • a) where the vector-charge 

- 2E a = (0, -2e ) (27) 

is the background charge in the Coulomb gas . The value of E is chosen in such a 
way so that only the B and C bonds acquire a phase (±7reo) when crossing the seam. 
As a reminder, we note that the BC loops in the three-colouring model were chosen 
earlier to represent the loops of the FPL model; this choice is of course arbitrary, 
and any one of the three possible choices of colour pairs would give the same results. 
Summing over the two possible ways of colouring a BC loop, we find for the fugacity 
of loops winding around the cylinder: 

n = 2 cos(7reo) , (28) 
which when compared to equation ( p6|) gives the relation 

7re = 6A . (29) 
The effective field theory in the presence of the seam can be written as [p6[ : 



/ = I d*rwg(\Vh 1 \ 2 + \Vh 2 \')+2itie (h 2 (L,oc)-(h 2 (L,-oc)).(30) 
This modification of the Gaussian field theory results in a shift of the conformal charge, 



e 2 



c=2-6^, (31) 
9 

and it also has an effect on the scaling dimensions of operators. As shown by 
Dotsenko and Fateev p3| , the correlation functions of the modified Coulomb gas 
can be expressed in terms of correlation functions in the Gaussian model (^|). The 
non-zero two-point functions in the modified theory are those of vertex operators 
having opposite electro-magnetic charge in combination with a floating charge 2E 
that precisely cancels the background charge —2E . The floating charge may combine 
with the negative charge. An operator whose total electro-magnetic vector-charge is 
(G, b) then has a dimension x(G, b) given by (ll| 

x(G,b) = 1g-(G-2 j B ) + ||6| 2 . (32) 

f Here we have assumed that the staggered chirality is an integrably marginal operator This, 
as will be shown later, is confirmed by Baxter's solution ]l4| . 
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Instead, the floating charge may split up so that both the positive and negative charge 
get increase by E . The dimension of such an operator is given by x(G + E ,b). This 
is the complete spectrum of scaling dimensions in the FPL model, where G G 7Z*, 
b G 1Z. The background charge is related to the fugacity by equation (|2|) , while the 
relation between the coupling g and the fugacity n remains unknown. We turn to this 
problem next. 

In the n = 2 FPL model the staggered chirality (or equivalently, the locking 
potential) is marginal. It follows from the exact solution obtained by Baxter [|l4| that 
it remains marginal as long as exp(i6A) lies on the unit circle, i.e. A is real. Therefore, 
the chirality is marginal in the entire regime < n < 2, and one expects it to change 
the value of the rcnormalizcd coupling constant g continuously with n. The dimension 
of the staggered chirality, for n < 2, is governed by the second smallest vector in 1Z* 
parallel to Eq. Using equation ( |32| ) with G x = (0,2), the marginality of x(r) (i.e., 
x(G x ) — 2) gives the relation between g and e§: 

1^-2. (33, 

where, from equation ( pi] ) we can read of the dependence of the background charge 
on the fugacity, 

eo = 1 — g = — arccos ( — ) . (34) 
7T V 2 / 

This agrees with the relation obtained numerically by Blote and Nienhuis |7j , and the 
Bethe Ansatz result of Batchelor et al. 

Equipped with equations (^), (H^), and (]34|) we can now calculate the watermelon 
dimensions x m for the n < 2 FPL model. In the Coulomb gas representation the 
watermelon scaling dimensions are given by x(Eo,b m ), where the magnetic charge is 
given by equation ( ft2|) for m even, and by equation (|2^) for m odd. The electric 
charge Eq is due to the electric- type operators exp(i27rii?o ■ h) that must be inserted 
at the endpoints of the watermelon configuration in order to correct for the spurious 
phase factors exp(±i6A) that arise due to the winding of the loop segments around 
the endpoints . We conclude that the watermelon scaling dimensions are given by 

/ft r % 5,2 (1 ~ af 

X2k = X{E ,b 2 k) = 7T « 7: 

2 2g 

x 2fc _! = x(E ,b 2 k-i) - f (fc 2 - fc + 1) - (1 2g g) , (35) 

which was also found by Batchelor et al, (see equations (16) and (17) in Q) and it 
generalizes the n — 2 result (^4j). The dimension x 2 was also calculated numerically 
by Blote and Nienhuis for different values of n, and their results agree very well 
with the exact values. 

The temperature dimension xt, for different values of n, was also determined 
numerically by Blote and Nienhuis, but this dimension does not appear in the Bethe 
Ansatz solution of Batchelor et al.. The dimension xt is related to a defect in the 
FPL model associated with an uncovered vertex. Unlike the case of the watermelon 
dimensions there are no loop segments associated with this defect and there is 
consequently no need for an electric type operator to correct for the winding of the 
loop segments around the endpoints. Therefore, we have 

3 

x t = x(0, bt) = -g , (36) 
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where the magnetic vector-charge b t = 3 A ((|T^)). The numerical results in Q are in 
very good agreement with our exact result. 

Once again the equations (pq) and (pq) are special cases of a more general spectrum 



x iJk = | (f + k * + jk) - (1 - S SJk ) (37) 



f L = — In A ~ /oc + 7771 . (39) 



for the charge bj t k = j(A — B) + k(A — C). These, like equation fl25|), govern the 
correlations between mixed defects defined by empty and occupied defect edges and 
the requirement that the occupied edges of both defects are mutually connected. Only 
for defects with topological charges purely in the A direction is the exponent unaffected 
by the background charge. 

Another way of viewing the effect that the background charge has on the exponents 
is the following. Consider the transfer matrix of the FPL model on a cylinder of 
circumference L. In terms of the eigenvalues of the transfer matrix (Ao > Ai > A2 > 
. . .) the exponents are defined by the gaps between the largest eigenvalue and the 
smaller ones |2q| . 

The free energy of the FPL model is given by the logarithm of the largest eigenvalue, 
and its finite size scaling is given by 

Z lnA ~ /oo + — 

where foo is the bulk free energy, and c the central charge. The background charge 
modifies the ground state of the system through the central charge, equation (|3l|), 
and therefore the various exponents, corresponding to excited states with respect to 
this shifted ground state, will also change. In the case of the watermelon dimensions 
there are extra lines running along the cylinder. This means that for this excited 
state we should "turn off" the seam, otherwise one would count the winding number 
of the loops around the cylinder which are not permitted, since they would intersect 
the lines along the cylinder. This state will therefore not have the correction from 
the background charge. As it is measured with respect to the new ground state via 
equation (|3g|) , the watermelon dimension x m will not only differ from the n = 2 value 
(x(b m ), equation (p4|)) due to the different value of g, but also by a shift, which is 
simply 1/12 of the shift of the central charge in equation (fjjj). Therefore, the value of 
the watermelon dimensions is 

Xra = X(b m ) - 6(1 ~ (40) 

IZg 

which coincides with equation j35|). In the state corresponding to the thermal 
excitation loops going around the cylinder are permitted and we need the seam to 
count those contributions correctly. This state therefore gets shifted in the same way 
as the ground state. The thermal exponent therefore only differs from the n = 2 
exponent in the value of g, and is given by equation (|36|). 



5. Geometrical Properties of Loops 



In this section we consider the geometrical properties of loops in the FPL model. We 
calculate the fractal dimension of loops and the loop length distribution function. 



15 



The length of a loop s, and its radius R are related by 

s ~ R D > , (41) 

where Df is the fractal dimension of the loop. The loop radius is defined as the radius 
of the smallest circle that contains the loop. 

The distribution of loop lengths P{s) measures the probability that a loop, in a 
fully packed configuration, passing through a chosen vertex has length s. It is given 
by a power law: 

P(s) ~ s-^- 1 ) . (42) 

The geometrical exponents Df and r can be related by a scaling argument to the 
exponent x 2 |27], ^ 

D f = 2 - x 2 , t 1 = — ?— . (43) 

Using the calculated value of x 2 , equation (|35|), we find: 

Dt = 1 + T 9 > r " 1 = iW (44) 

In the limit n — > (g = 1/2) we find Df = 2, which is what we expect for 
Hamiltonian walks. For n — 2 we find Df — 3/2 which is the fractal dimension of 
equal height (contour) loops on a random Gaussian surface |24| . Indeed, for n — 2 the 
background charge is equal to zero and the free energy / defines a random Gaussian 
surface; the loops in the FPL model are contour loops of this random surface. 

In the Kagome Potts model representation of the n = 2 FPL model the exponent 
t — 1 has been determined numerically to be 1.34 ± 0.02 |28||, in good agreement with 
the exact result r — 1 = 4/3 which follows from equation (f44j) for g = 1. 



6. Summary and remarks 



We have calculated the conformal charge and the exact exponents in the fully packed 
loop model on the honeycomb lattice. To this end we proposed a simple conformal 
field theory for the scaling limit of this model, and it is given by the vacuum phase of 
a two dimensional Coulomb gas with an added background charge. The magnetic and 
electric charges of the Coulomb gas were found to be vectors in the triangular lattice 
1Z (the so-called "repeat" lattice) and its dual TZ*. These charges give the complete 
operator spectrum of the fully packed loop model of which the watermelon dimensions 
calculated by Batchclor et al. are a subset. The exact value of the temperature 
dimension found here and not in the Bethe Ansatz solution are in agreement with the 
numerical results of Blote and Nienhuis . 

Coulomb gas methods with vector-charges have been used previously by Fateev 
and Zamolodchikov [^9| to calculate correlation functions in the Z3 models. Their 
work was extended by Pasquier ]30| who considered the continuum limit of lattice 
models with quantum group symmetries. The FPL model is most likely related to the 
models of Fateev and Zamolodchikov since in the n — > 2 limit the FPL model has an 
enlarged chiral symmetry, given by the su(3)k=i Kac-Moody algebra This is also 
true of the Z3 models. Moreover, for n < 2 we introduced a term proportional to the 
staggered chirality into the free-energy (action) of the FPL model. The presence of 
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this term breaks the full permutation symmetry of the three-colouring model down to 
Z3; only cyclic permutations of the colours leave the staggered chirality unchanged. 

Finally, is interesting to note that in the folding model the microscopic heights 
z(r) specify the positions of the vertices, r, of the triangular lattice in the folded state. 
The coarse grained (entropic in origin) free energy, corresponding to the different 
ways of folding the triangular lattice, is Gaussian, equation (||). In the theoretical 
considerations of tethered membranes this is usually the starting assumption which is 
justified by the results of numerical simulations |3^]. We see that in this simple folding 
model the gradient squared form of the entropic contribution to the free energy of 
folding, is closely related to the fact that the FPL model maps to the vacuum phase 
of a two-dimensional Coulomb gas. 

The ideal states we had identified in the three colouring model are the states that 
are entropically selected. They are flat, h = const, and in the folding model they 
map to states in which the whole triangular lattice has been folded down to a single 
triangle. Introducing an energy penalty associated with folding might stabilize a flat 
state of the membrane at low temperatures, while the entropically selected folded 
states would necessarily win at high temperatures. Therefore, we might expect a 
folding transition to occur at some intermediate temperature. Such a transition was 



found in the numerical transfer matrix calculation of Di Francesco and Guitter 16 . It 
would be interesting to study this transition using the Coulomb gas methods developed 
here. It is our hope that this will allow us to calculate properties of this intriguing 
transition exactly. 
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